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Abstract 

We investigate quintom cosmology in FRW universes using isomorphic models consist- 
ing of three coupled oscillators, one of which carries negative kinetic energy. In particu- 
lar, we examine the cosmological paradigms of minimally-coupled massless quintom, of two 
conformally-coupled massive scalars and of conformally-coupled massive quintom, and we 
obtain their qualitative characteristics as well as their quantitative asymptotic behavior. 
For open or flat geometries, we find that, independently of the specific initial conditions, the 
universe is always led to an eternal expansion. 
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1 Introduction 



The type la supernova observations suggests that the universe is dominated by dark en- 
ergy with negative pressure, which provides the dynamical mechanism for its accelerating 
expansion [2 [3] [3]. The strength of this acceleration is presently a matter of debate, mainly 
because it depends on the theoretical model implied when interpreting the data. 

The most obvious theoretical candidate for dark energy is the cosmological constant A (or 
vacuum energy) which has the equation of state w = — 1. However, as it is well known, 

there are two difficulties arising from the cosmological constant scenario, namely the two 
famous cosmological constant problems — the "fine-tuning" and the "cosmic coincidence" 
one [BJ. An alternative proposal is the concept of dynamical dark energy. Such a scenario is 
often realized by some scalar field mechanism and suggests that the energy form with negative 
pressure is provided by a scalar field evolving under a properly constructed potential. So far, 
a large class of scalar- field dark energy models have been studied, including quintessence [7J, 
K-essence [8], tachyon [9], phantom [10], ghost condensate quintom |12] , and so forth. 
It should be noted that the usual viewpoint regards scalar-field dark energy models as an 
effective description of an underlying theory of dark energy. In addition, other proposals on 
dark energy include interacting dark energy models [13] . braneworld models [14] . Chaplygin 
gas models [15], holographic dark energy [16], bulk holographic dark energy [17] and many 
others. In this context, scalar fields, which may come in different forms and with a variety 
of possible self-interaction potentials, constitute the dominant (sole) form of matter in the 
field equations of the gravitational theory. 

In the present work we consider scalar fields as the ones responsible for a dark energy uni- 
verse behaving as quintom, that is obtained from an interplay of phantom and quintessence 
models. Neither of these two models alone can fulfill the transition from w > — 1 to ra < - 1 
and vice versa. Furthermore, although in k-essence [3] one can have both w > —1 and 
w < —1, it has been lately shown in [181 1191 that the corresponding crossing is very unlikely 
to be realized during the evolution. However, one can show [12U20] that considering the com- 
bination of quintessence and phantom in a qualitatively new model, the — 1-transition can 
be fulfilled, as can be clearly seen in [12]. The quintom scenario of dark energy is designed 
to understand the nature of dark energy with w across -1. The quintom models of dark 
energy differ from the quintessence, phantom and k-essence and so on in the determination 
of the cosmological evolution and the fate of the universe. 

Under the assumption of a quintom field with negative kinetic energy, the demand for 
a final dominance of phantom universe leads naturally to the consideration of a minimally- 
coupled quintom field, too. In the non-spatially flat universe, the problem of quintom stabil- 
ity has motivated the formulation of a toy-model [21] consisting of three coupled oscillators, 
one with negative and the others with positive-definite kinetic energy, with the first oscil- 
lator mimicking the gravitational field while the other two mimic the quintom (phantom 
and quintessence) field. In section [2] we summarize the equations for the coupled-oscillator 
models. In particular in subsection 12.11 we formulate the dynamics of the massless quintom 
scenario, while in 12. 21 we construct an oscillator-model with positive energy which describes 
exactly a conformally-coupled quintom in rescaled variables and conformal time. In section 
[3] we investigate the characteristics of the isomorphic model of the massless quintom and in 
section [4] of that corresponding to the conformally-coupled quintom. Finally, in section [5] we 
summarize the obtained results and we discuss their physical implications. 

2 Coupled oscillators as cosmological models 

We know that toy-models, in which coupled oscillators are utilized, are present in recent 
cosmological studies. In this section we rewrite the field equations relating the quintom 
model and scalar field cosmology and we show that these equations can be reduced to the 
equations of a model with three coupled oscillators. We perform our investigation in the 
framework of a non-spatially-flat Friedmann-Robertson- Walker (FRW) universe which is 
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described by the line element 



ds 2 = -dt 2 + a 2 (t) { dr - + r 2 dfl 2 I . (1) 
1 — kr 2 



2.1 Quintom cosmology 

In order to describe the quintom field we use, as usual, two scalars: <j> and a. For simplicity 
we assume that they constitute the only form of matter present. The energy density and 
pressure of this minimally-coupled quintom field in the metric © are respectively: 



1 12 1 . 2 

— q> H — a 
T 2 



P = 



2 . 1.2 
> +-C 



V(<t>,a) 



(2) 



(3) 



which correspond to the normal definitions for minimally-coupled scalars but with sign- 
inversion of the kinetic energy of <f>. In a homogeneous and isotropic universe <f>(t) and a(t) 
depend only on the comoving time t and interact through the potential V{4>, a). V(<j>,a) 
should not be more negative than a norm because, due to the final dominance of phantom 
universe, the energy density p should be non-negative and the Hubble parameter H should 
be real. Note that a minimally-coupled scalar with positive kinetic energy density + h4> 2 
cannot exhibit P < —p in Einstein gravity [22] . 

For a quintom-universe described by metric ©, with the 4> an d o fields as material 
sources, the equations of motion are the following: 



-j> 2 + 2V(<P,a)] 



H + H 2 = ^ [~£r 2 + 2 + V(cP,*) 
dV 



(4) 
(5) 
(6) 



dV 

a + 3H& + ^~ 
aa 



0, 



(7) 



where k = 8ttG and dots denote differentiation with respect to the comoving time t. It is 
remarkable that only three of the equations (©-(El are independent. Indeed, when <j> 7^ or 
it / one can derive one of the Klein-Gordon equations (O,© using the other and either 
(jl]),© or the conservation equation p+3H(P+p) = satisfied by the quintom. Furthermore, 
a significant feature of equations (©-©i which determine the dynamical evolution of the 
universe at hand, is that ©, © and a combination of © and ©, can be derived from the 
Lagrangian: 



Lo = Ka 3 (p — P) = 3ad 2 — 3afc + (ca 3 



?-\cr 2 + V{4>,o) 



or from the Hamiltonian: 
Ho = 



3ci 



H — — + — 



(8) 



(9) 



The aforementioned simple quintom cosmological model comprises of three degrees of 
freedom, one of which carries negative kinetic energy. Therefore, we can construct a toy- 
model that consists of three coupled oscillators, one with negative-definite and two with 
positive-definite kinetic energy. This toy-model, although not exact (i.e re-producing the 
true system), it presents the same qualitative behavior [23] and thus it can mimic the real 
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system described by © or ©, in the case of k = 0. In the simple case of a quadratic 
potential V = m 2 (4> 2 + a 2 )/2, such a model can be formulated using the Lagrangian: 



•2 -2 22 2 2 222 \222 

2/ + £ m^a: _ m y y _ m z z _ fj, x y _ X x z 



•2 -2 „,2 2 „, 2 „ Jl , 2„2„ 2 

SL + i. 

2 2 2 
or the associated Hamiltonian: 

■2 -2 -2 22 2 2 222 \222 

x v z m^x m v y m z z a x y X x z , 

ll = Y~-2 + Y + ^ + ^r + — + tL ^ L + ^- 

As far as the total energy of the system remains constant, the energy of the phantom oscil- 
lator could decrease arbitrarily, while the energy of the other two oscillators could increase 
infinitely. Hence, there is not a stable ground state for the system [231 121j . 
We can now derive the Euler-Lagrange equations from (|10|l : 

x + (m 2 + i-i 2 y 2 + X 2 z 2 )x = (12) 

y - (m% + n 2 x 2 )y = (13) 

z + (m 2 z + X 2 x 2 )z = 0. (14) 

In the following we assume that m x — m y — m z = and fi 2 — 1 and A 2 = 1, so that 
equations (|12[) - (|14p are reduced to: 

x = -(y 2 + z 2 )x (15) 

y = x 2 y (16) 

z — —x z. (17) 

The toy-model governed by the equations of motion (115[l - (|17[ l mimics the behavior of that 
of l|4" |l -([? |) . or in other words the massless quintom cosmological paradigm is qualitatively 
isomorphic to the constructed system of coupled oscillators. This property allows us to reveal 
the characteristics of the cosmological scenario by investigating the dynamical evolution of 
this oscillator system. 

2.2 Conformally-coupled scalar field cosmology 

There are many arguments supporting that in a curved space a scalar field couples non- 
minimally to the Ricci curvature. The explicit non-minimal coupling to the curvature intro- 
duces extra terms in the equations for the scalar fields and the scale factor, allowing for an 
accelerating expansion for the universe or even for a super-accelerating (H > 0) universe. 
Apart form the possibility of explaining the observed recent universe acceleration, there are 
additional reasons to consider such a model. Non-minimal coupling is introduced by quan- 
tum corrections to the action of a scalar field [24], it is necessary for the renormalization 
of the scalar field theory [25] . and it is even required at the classical level to preserve the 
Einstein equivalence principle or to avoid causal pathologies Thus, in this work we 

consider conformally-coupled massive scalar fields. 

We use an action including scalar fields with positive, non-minimally-coupled kinetic 
energy with the Ricci curvature [23] : 



j4 

a X\ 
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Y K - \ + ° 2 ) R ~ i<? ab V^V^ - \g ab V a aV b o - V(<t>, o) 



(18) 



where £ is a dimensionless coupling constant. We are interested in studying the specific case 
where f = 1/6, since such a choice is an infrared fixed point of the relevant renormalization 
group [27] . Assuming a potential form V = m 2 ((j> 2 + a 2 )/2, we derive the equations of motion 

Eni ES [321 US]: 



,2 k Km 2 

T -1 

a 2 6 



H + 2H 2 + ^-^U 2 +a 2 )=0 (19) 
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|(0 2 - a 2 ) + k,H{H - a&) - 3H 2 



|(^ + «r a )J+^(^ + 0=° ( 20 ) 

<f> + 3H<f>+ —4> + m 2 (j> = (21) 
6 

cf + 3if<j - -a - m 2 a = 0, (22) 
6 

along with the Hamiltonian constraint as the first fundamental FRW equation: 

" 2 +^ = fp- (23) 

The effective energy density and pressure of the scalar fields, which guarantee energy con- 
servation, are given by |34j : 

p = i (0 2 - a 2 ) + ?L (0 2 + a 2 ) + JL (0 2 + a 2 ) + itf0(#0 + 20) + ijfa(H<r - 2a) (24) 
and 

p = i (0 2 - . 2 ) - f (0 2 + - £ (0 2 + ^) - i [m* - -)+ 

+2(0 2 - a 2 ) + 200 - Igo + (2H + 3// 2 )(0 2 + a 2 )] . (25) 

Equations (I19|l - (|22| ) and expression (|24p for p are evidently complicated. However, we can 
reduce the problem to a system of three coupled oscillators with sharply defined energies in 
a fictitious Minkowski space, repeating the same steps as in the previous subsection. Indeed, 
inserting the auxiliary variables: 



introducing the rescaled (conformal) time r\ given as dt = adrj, and assuming a quadratic 
potential V = m 2 ((j> 2 + a 2 )/2, the equations of motion are transformed to [331 1361 13"7]: 

x" = (y 2 + z 2 - k)x (27) 

y" = ~x 2 y (28) 

z" = x 2 z, (29) 

where primes denote differentiation with respect to r\. Note that these equations can be 
obtained from the Lagrangian: 

L i = -\( x 'f + \(y'f - \( z ')' 2 - \ x ' 2 y 2 - \ x2z2 + \ kx ^ ( 30 ) 

or from the Hamiltonian: 

HI/ /\2 . 1/ /\2 1/ f\2 .122,122 lj 2 / 01 \ 
i = -^\ x ) + TjlW ) - ) + 2^ y + 2^ ^ ~ 2 ' *■ ^ 



Clearly, the system of equations ((27j - (J29J is isomorphic to that of (|19[) - (|22|) . We mention here 
that y and z are a mixture of the gravitational and scalar-field degrees of freedom, whereas 
x is associated solely with gravity. Furthermore, as it is required by the first relation of (|26|l , 
the restriction x > must be applied. 
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3 Dynamical behavior of massless-quintom oscillator- 
model 



In the previous section we formulated the use of oscillator-models in investigating the dy- 
namical behavior of cosmological paradigms. For a spatially flat universe under the massless 
phantom scenario, the associated system of two coupled oscillators was studied by Castagnino 
el al [37] (see also [23]). In this section we examine the evolution characteristics of the 
oscillator-model corresponding to the massless quintom, which consists of three degrees of 
freedom and was formulated in subsection l2.ll In this case the phase-space stability analysis 
leads to the following results: 

-Firstly, we extract the fixed points of the system (|15j> - (|17[l . defined as those points 
where the velocities of the oscillators are zero. It is easy to see that these are simply the 
loci (xo, 0, 0) and (0, yo, zo). In order to examine the stability of these fixed points, as usual 
we calculate the partial derivatives of the right hand sides of the system (|15|l - (|17| ) at these 
points, and we extract the eigenvalues of the corresponding matrix [3S]. Since at least one 
of the eigenvalues is always positive, we conclude that the fixed points are all unstable. In 
fact, a three-oscillator system with one of them having negative-definite kinetic energy, does 
not possess positions of equilibrium. 

-Due to these instabilities, apart from the fixed points, all the orbits in the phase-space 
go to infinity as t ^ oo. In particular, x(t) — > oo monotonically while y(t) and z(t) — » 
oscillating, a behavior which is independent of the choice of initial conditions. This can be 
easily verified by simple numerical investigations. 

-Cycles (periodic orbits) are not possible. However, chaotic dynamics may appear. This 
is a robust result for the case of the quadratic (or equivalently the Yang-Mills) potential, 
and arises from the corresponding extensive studies of the literature (see for example 39 ). 

-In general, as it has been shown in [231 137] . and taking into account the invariance under 
the transformation (x,y,z) — ► (— x,y, z) and (x,y,z) — > (x,—y,z) and (x,y,z) — > (x,y, — z), 
we can consider an asymptotic solution as: 

V(t) «*(t)«^ sin (~\ (32) 

x(t)&t 2 . (33) 
In this case the kinetic energy of the y and 2-oscillators (for large times) is 

K M = (£± (34) 

We can see that this expression oscillates with divergent amplitude, while the kinetic en- 
ergy of the :r-oscillator, = — ~ — 2t 2 — > — oo. This behavior corresponds to the 
instability described in [21] . 

Having examined the system characteristics using the auxiliary degrees of freedom (|26[l 
we can now transform back to the physical variables a, (j> and a, under the restriction x > 0, 
noting the necessary inversion between variables indicated in [23] . Doing so our results can 
be re-written as follows: 

-The fixed points of the cosmological model are just (do, 0,0) (we discard the family 
(0,j/o, Zo) since it corresponds to the non-physical case of a universe with zero scale factor). 
Thus, they correspond to Minkowski spaces with constant scale factor, i.e without expansion. 
The fact that they are unstable, i.e the absence of attractors, implies that an arbitrary small 
perturbation can lead to the aforementioned diverging orbits. 

-In particular, the previously analyzed behavior of x(t), y(t) and z(t) implies that the 
scale factor a goes to infinity, while the scalar fields are oscillating, in general with a varying 
period. Thus, we conclude that we always acquire an expanding universe. This is true 
even in the case where the system lies initially in one of the fixed points (constant scale 
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factor and zero derivative), since an arbitrary small perturbation is sufficient to lead it to 
the aforementioned expanding case. Finally, as was shown in (25], the asymptotic solutions 
(|32|) . (|33|) correspond to a matter-dominated universe with a(t) = aot 2 ^ 3 . That is, in this 
model there is no mechanism that can end the expansion, either by reversing it to contraction 
or by stabilizing the universe to a steady-state type. Thus, increasing dilution and the 
"thermal death" of the universe are inevitable. 

-Periodic orbits do not exist. That is, a massless quintom cannot drive cyclic universes 
[40] . However, since chaotic behavior is possible, we conclude that we can obtain chaotic 
cosmological evolution. Indeed, FRW cosmologies are known to present chaotic behavior 

an- 



4 Dynamical behavior of conformally-coupled quin- 
tom oscillator-model 

Let us now examine the dynamical characteristics of a massive conformally-coupled quin- 
tom, extending the associated oscillator-model formulated in subsection 12.21 (see [42] for the 
corresponding problem for a phantom field). The Klein-Gordon equations of the scalar fields 
are: 

4> + 3H<p-m 2 (j>-tB4> = (35) 

a + 3.Hct + raa + £Ra = 0. (36) 

In terms of the auxiliary variables x, y and z defined in (|26p and the conformal time dt — adr], 
these field equations can be derived from the Lagrangian: 

12 12 12 2 2 2 2 i 

x y z x y x z 1, 2 

or from the Hamiltonian: 

12 12 12 _2 2 _2 2 i 

x y z x y x z 1,2 /„„n 

H 2 = h — + h -kx 2 . (38) 

2 2 2 2 2 2 y J 

Note that the Lagrangian (|37[) is equivalent to Li of equation (|30[) (apart from an overall 
sign) provided that the conformally-coupled scalar is turned into a phantom field. The 
dynamical system in this case is: 

x" = {y 2 + z 2 - k)x (39) 

y" = x 2 y (40) 
z" = x 2 z. (41) 

In order to perform a stability analysis, we have to distinguish between the various fc-cases. 

-Firstly, in the case of a flat universe (A; = 0) we can see that the fixed points are just 
(xo, 0, 0), while for k = 1 they are (xq, ±f, ±\/l — v 2 ). Both these loci correspond to empty 
Minkowski spaces. However, for the case of a closed universe (k — —1) there are not fixed 
points, apart from the trivial case of a zero-scale-factor universe. 

-It is straightforward to see that these fixed points are unstable, and that all the (non- 
stationary) orbits in the phase-space go to infinity as t — > oo .37]. In particular, numerical 
investigation shows that for k = 0, 1, x(t) — > oo monotonically while y(t) and z(t) — > 
oscillating. In fact, the oscillatory nature of the solutions for k — is a general feature of 
non-minimally coupled scalar fields with £ > 0, even in the case where £ 7^ 1/6 [30]. Finally, 
in the case k = — 1 there is not a qualitatively general asymptotic behavior and the system 
evolution can be more complicated. In this case, as is confirmed by numerical integration, 
the scale factor can sometimes decrease. 

-Thus, in the cases of open and flat universes, all the solutions that are not stationary 
represent universes expanding to infinity. However, due to the instabilities, we conclude that 
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even an arbitrary small perturbation can bring the universe out of stationarity, and lead it 
to an everlasting expansion. Therefore, independently of the initial conditions, we always 
obtain an expanding universe. The expansion cannot be reversed or end, and a complete 
dilution is inevitable. Finally, note that contrary to the case of the previous section, there 
is not a specific quantitatively asymptotic behavior. Thus, the Hubble parameter can be 
constant, or increasing, leading to an accelerating or even super-accelerating universe. On 
the other hand, for the case of a closed universe, the system evolution is more complex and 
difficult to be outlined. 

-Closed orbits in phase-space do not exist, that is a massive conformally-coupled quintom 
cannot bring about a cyclic universe |40j . However, due to the quadratic nature of the 
potential, chaotic behavior is possible, which, as we have already mentioned, is expected in 
FRW cosmologies [41 J . Finally, chaoticity is more easily obtained in the case of a closed 
universe (k = — 1), where the dynamics of the system can be more complicated. 

5 Conclusions 

In this work we investigate the evolution characteristics of quintom universes, using oscillator- 
models as an isomorphic description. Indeed, one can construct such models which present 
all the relevant information, examine their dynamical behavior under stability analysis, and 
finally transform the results back to the cosmological picture. In particular, in [37] the au- 
thors have constructed a toy-model with two degrees of freedom in order to isomorphically 
describe the cosmological paradigm of a minimally-coupled massless phantom field and a 
massless graviton in physical time. In such a case the perturbations of the model are un- 
stable and this feature led the authors of |21) to extend it by the insertion of a negative 
mass term to the potential, necessary for the stabilization of perturbations (the model of 
Castagnino el al. of 37 could be recovered by forcing the phantom to be massless). 

The dynamical system (|15[l - (|17[l of the present work constitutes a toy-model for a mini- 
mally coupled massless quintom field and a massless graviton in physical time, since in this 
case the isomorphic description requires three degrees of freedom. Expressing the results in 
the physical framework we find that the fixed points of the system are just the Minkowski 
spaces, but are unstable under, even arbitrarily small, perturbations. Moreover, there are 
no attractor points and all the phase-space orbits go to infinity with increasing time. In par- 
ticular, in this asymptotic case we obtain a matter-dominated universe with a{t) = dot 2 ' 3 . 
Thus, we conclude that, independently of the initial conditions, the universe is always led 
to an eternal expansion. 

The dynamical system (|2T[> - 129[) describes two conformally-coupled massive scalar fields. 
For flat and open geometries, all the fixed points are empty Minkowski spaces and they are 
unstable. With increasing time the scalar fields go to zero in an oscillatory way. The obtained 
evolution reveals that the matter content dilutes progressively up to its complete evanescence 
while the universe expands. Such a behavior is consistent with the phenomenology of non- 
minimally-coupled scalar fields. In addition, periodic orbits are not present, and thus this 
cosmological paradigm cannot drive cyclic universes. However, chaotic dynamics may arise, 
a feature already known for FRW cosmologies. On the other hand, for a closed universe, 
the dynamics of the system is more complex, and it is hard to determine even a general 
qualitative behavior. 

The dynamical system (|39|l - (|4ip corresponds to a conformally-coupled massive quintom 
field. Investigating the behavior of its phase-space we see that for a flat or open universe the 
fixed points represent Minkowski spaces, which are unstable under perturbations. Thus the 
universe is expanding to infinity. However, this scenario allows for a general (unspecified) 
Hubble parameter, i.e for either accelerating or even super-accelerating universe. Closed 
orbits, corresponding to cyclic behavior, do not appear but chaoticity does. Finally, for a 
closed universe the system evolution can be more complicated. 

In conclusion, we observe that we can extract qualitative as well as (asymptotically) 



quantitative characteristics of various quintom paradigms, by investigating the corresponding 
isomorphic coupled-oscillator models. The stability analysis and the obtained asymptotic 
behaviors show that the quintom scenario is consistent with cosmological observations. 
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